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1. Introduction 

The family of carbon-based materials, including two distinct crystal forms: the isotropic 
diamond and anisotropic graphite, together with other allotropes such as fullerenes and 
carbon nanotubes, has been recently enlarged to include graphene, a unique material 
which consists of a two-dimensional (2D) lattice of carbon atoms with a honeycomb 
symmetry. Graphene, the closest archetype of a two- dimensional crystal, while being 
very promising for future applications, represents also an extremely interesting system 
from the viewpoint of fundamental physics [H |2]. Distinctly, the electronic states of 
graphene do not obey the conventional laws of Schrodinger's quantum mechanics but 
are rather governed by an equation which is equivalent to the Dirac equation for massless 
fermions [21 H]. This very peculiar, but at the same time simple, electronic band structure 
of graphene has been investigated theoretically for more then sixty years [5]. The "rise" 
of graphene dates however only from 2004 when it was effectively introduced into the 
the laboratory environment, following the clear identification of 10 /im sized graphene 
flakes deposited on Si/Si02 substrates [6] and of multi-layer epitaxial graphene thermally 
decomposed from silicon carbide [71 [8] . More recently, the growth of large area graphene 
samples on metallic surfaces (Ru, Ni or Cu) [HI [ini [HI [121 [13 [E] , and their subsequent 
transfer on arbitrary substrates has also substantially progressed [HI HI], what may be 
a promising option for easy preparation of macroscopic-size crystallites of graphene. 

The intense, over last years, investigations of different graphene systems have 
expanded today into a broad research field of qualitatively new two-dimensional 
systems [15]. The challenging expectation are that graphene will introduce quantum 
electrodynamics into solid state laboratories [TH [T71 HHl [19] and that it will constitute the 
basis for future electronics plausible avoiding the known limitations (at the nano-scale) 
of the CMOS technology. As for today, the basic properties of "new" 2D allotropes of 
carbon, including graphene [31 S], graphene bilayer [201 121], multi-layer graphene [6[ 122], 
graphene on a silicon carbide substrate [7], are already relatively well-known and the 
basis of the graphene physics is pretty well established. 

Despite relevant advances in various experimental techniques, such as angular 
resolved photoemission spectroscopy (ARPES) [231 [211 [251 [26] or scanning tunnelling 
spectroscopy (STS) [23 [2H1 [29], a large part of our knowledge on electronic properties 
of graphene has been deduced from conventional, electric transport and optical 
experiments, in particular under an applied magnetic field. The physics of the integer 
quantum Hall effects in single-layer and bilayer graphene [3l [H [20] or of the minimum 
conductivity of graphene [30l [31] is nowadays underpinned by a solid basis of published 
transport experiments. Optical investigations of graphene systems are equally fertile. 
They offer a rather direct verification of electronic band structure in these materials 
and/or help to avoid the advanced sample processing of often invasive character [32] . 

A number of papers which review the graphene-related physics has already appeared 
within last few years [33l [H [HH [2l [35l [36l [15], but the progress in this research area 
is extremely fast. This, in particular, concerns the developments related to optical 
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spectroscopy which was perhaps not very popular at the birth of the research on 
graphene but it is now more and more frequently used, as testified by a large number of 
recent papers cited here. As any review article also this one includes some unavoidable 
repetitions of information which can be found elsewhere. However, our intention is 
to discuss the properties of different graphene systems from the viewpoint of optical 
spectroscopy methods, including the most recent developments in this field. Whereas 
the previous review articles were mostly focused on the properties of a single graphene 
monolayer, here we discuss the Dirac-like electronic states (massless and massive) which 
in fact are representative of any two-dimensional allotropes of carbon or even of bulk 
graphite. 

Indeed, the methods of optical spectroscopy have been recently applied to a 
large class of graphene-based systems. A number of various optical and magneto- 
optical measurements have been performed on mono-, bi- and multi-layer Bernal-stacked 
graphene (exfoliated from graphite and placed on Si/Si02), on multilayer epitaxial 
graphene thermally decomposed from C-terminated surface of SiC, as well as on bulk 
graphite. The published papers provide a mosaic of information as they usually address 
a particular problem or individual sample characteristics. Investigations of the electronic 
band structure remains, however, the key element of the majority of these papers and 
also constitutes the bulk of our review. Other characteristic problems invoked with 
optical methods concern the efficiency or mechanisms of carrier scattering. The effects 
of disorder or interactions are also here discussed. 

The paper is organized as follows. The next section has still some introductory 
character in which we stress the relevance of the optical spectroscopy in the "rise" 
of graphene, i.e., in the identification of the graphene fiakes on Si/Si02 substrates. 
The following sections 3-6 are systematically devoted to: 3) monolayer graphene and 
multilayer epitaxial graphene (on C-face of SiC) which both exhibit the simple Dirac-like 
electronic spectrum, 4) bilayer graphene with the unique massive Dirac-like fermions 5) 
bulk graphite, with more complex electronic bands but also displaying the Dirac-like 
dispersion relations, and 6) A^-layer (A^ > 2) segments of graphite. In each of these 
section we first present the theoretically expected band structure and optical response 
of the investigated system, then we confront it with the available experimental data 
and discuss the implications. The particular attention is focused on magneto-optical 
investigations, which we consider to be particularly valuable in the studies of electronic 
properties of the investigated systems. Conclusions and possible perspectives related to 
optical studies of graphene systems are presented in Section 7. 

2. How to see graphene 

The identification of the "exfoliated graphene" on Si/Si02 substrates by the group at 
the Manchester University [6] is one of the unquestionable milestone in the development 
of the graphene oriented research. Initially, however, it was perhaps not so clear why 
a few- or even single-graphene sheets can be seen just through an optical microscope. 
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Figure 1. Graphene crystallites on the Si/Si02 substrate with the oxide thickness 
of 300 nni imaged in white light (a), green light (b) and another graphene sample on 
the Si substrate with 200-nm-thick Si02 layer imaged in white light (c). The smaller 
figures in the upper and lower lines represent images taken at various wavelengths for 
the thickness of oxide layer 300 and 200 nm, respectively. The theoretically calculated 
contrast (for graphene monolayer) as a function of the wavelength and the thickness 
of the oxide layer is shown in part (d). In a relatively good agreement with these 
calculations, graphene crystallites are best visible in images taken at wavelengths 
slightly below 600 nm and around 400 nm, for the thickness of dioxide layer 300 and 
200 nm, respectively. Reprinted with permission from [37] . Copyright 2007, American 
Institute of Physics. 



The explanation of this today widely accepted and frequently used experience is not 
trivial and invokes the importance of the spectroscopic analysis already at the birth of 
the graphene physics. As more recently revealed in details [371 EH EH HOl SI] , the choice 
of the appropriate thickness of the Si02 layer plays the essential role in the visualization 
of the graphene layers. As shown in Fig. [H the contrast which is necessary to "see" the 
individual graphene layers depends on the wavelength of light and the actual thickness of 
the Si02 layer. This is accounted by the analysis of multireflections from the sandwiched 
graphene/Si02/Si structure within the framework of standard Fresnel equations (despite 
the fact that we are dealing with the material at the atomic resolution). The substrates 
used by Manchester group were characterized by 300-nm-thick oxide layers, for which 
the optimum contrast to see the graphene flakes is ~10 % and falls in the middle of the 
visible range. This greatly helped to visualize the graphene flakes. 

3. Graphene 

3.1. Electronic band structure 

The first calculation of electronic states in a 2D lattice of carbon atoms with a 
honeycomb symmetry, see Fig. [21 dates back to 1947, when Wallace [5] used graphene 
as a starting element for description of bands in bulk graphite. Taking into account the 
strong hybridization of 2sp^ orbitals in the graphene plane, Wallace considered just the 
remaining p orbital (oriented perpendicular to the crystal plane) to be responsible for 
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Figure 2. Segment of the graphene crystal: Carbon atoms are arranged into a 2D 
lattice with a characteristic honeycomb symmetry and lattice constant ap = 0.246 nm 
[42] , which is by a factor of larger than inter-atomic distance. 



the electronic band structure in the vicinity of the Fermi level and suggested a standard 
tight-binding approach. Considering only the nearest-neighbour hopping parameter 70 
(see Fig. [2]) , one easily finds a pair of 7r-bands [2] : 

E^*{k) = -EM = 



7o^ 



^1 + 4 cos^ + 4 cos cos , (1) 

which distinctly cross (touch) at two inequivalent K and K' points of the Brillouin zone, 
see Fig. [31 The strength of the nearest-neighbour hopping is 70 ~ 3.2 eV [13] and the 
lattice constant Qq = 0.246 nm |12] is by a factor of y/S larger than the distance between 
the nearest carbon atoms. 

In pristine graphene, the Fermi level lies just at the touching (crossing) point of vr 
and TT* bands and graphene has a character of zero-band-gap semiconductor (semimetal). 
Close to a given crossing (touching) point, known also as the Dirac or charge neutrality 
point, the electronic bands are nearly linear and practically rotationally symmetric. In 
other words, the carrier dispersion relations take a simple form: i?^* = — i?7r ~ "^F^lkl, 
where the momentum k is measured with respect to K {K') point. The parameter vp, 
having dimension of a velocity, is directly related to the coupling strength (hopping 
integral) between the nearest carbon atoms: Vp = y/Sao'jo/ (2h) [2]. 

The linearity of bands in graphene (in the vicinity of the K and K' points) implies 
that charge carriers in this material behave as relativistic particles with zero rest mass 
and constant velocity vp ~ 10^ m.s~^. They are often referred to as massless Dirac 
fermions, and with a good precision, their behaviour is described by the effective 
Hamiltonian [21 [S] : 

U f Px-iPy\ f 71^ \ , . 

H = vf[ n = ] =Vf(t-P, 2 

which is equivalent to the Hamiltonian in the Weyl equation for real relativistic particles 
with zero rest mass (originally for neutrinos) derived from the Dirac equation. Due to 
this formal similarity, a direct link between quantum electrodynamics and the physics 
of graphene is established. 
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Figure 3. Schematic view of the graphene band structure with the characteristic 
Dirac cones in the vicinity of K and K' points. 

The characteristic hnear dispersion relations of electronic states makes graphene 
very distinct among other 2D systems (such as quantum wells or heteroj unctions) which 
have been widely investigated in condensed-matter physics for last thirty years. For 
instance, the density of states in graphene is not constant as for conventional massive 
particles, D = g^gsm/ (Q.nh'^^, but rises linearly with the energy distance e from the 
Dirac point: D^e) = gsgy\e\ / {27iVph'^) . Here, gs = gv = stand for the spin and valley 
degeneracies, respectively. 

Nevertheless, it should be noticed that the relativistic-like image of electronic 
states in graphene given by Hamiltonian (jj]) remains an approximate model. This 
simple approximation is well, and perhaps even surprisingly well fulfilled in case of 
electronic states in the vicinity of the Dirac point. Those states indeed show almost 
perfect character of massless Dirac fermions when probed with different experiments. 
However, as can be seen from Fig. [3], the deviations from this relativistic model become 
obviously important for states far away from the Dirac point, even if we consider only 
the nearest neighbours in the tight-binding calculation. Other complications may arise 
when including the hopping integrals between further located neighbours. For example, 
when taking into account the non-zero values of 7o'^ hopping integrals |l5l HB] (see 
Fig. |2]), the nonlinearity is enhanced and Dirac cones become asymmetric with respect 
to the charge neutrality point [171 SHI HH]- 
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Figure 4. (a) Photograph of a 50-mni aperture partiahy covered by graphene and its 
bilayer, with scheraatically depicted optical absorption, (b) Transmittance spectrum 
of single-layer graphene in comparison with theoretical transmission for ideal massless 
Dirac fermions and real graphene. Inset shows relative absorption as a function of 
number of layers in graphite stacks. Reprinted from |49] with permission from AAAS. 

3.2. Optical conductivity at B = 

Interaction of massless Dirac fermions in graphene with hght, expressed usually in terms 
of dynamical (optical) conductivity, has been theoretically treated at different levels of 
approximation [50l[5lllll|52l|53l[Ml[55l|56lll6l[53 Within the simplest approach 

of electronic states described by Hamiltonian ([2]), only vertical in fc-space transitions 
across the Dirac point are optically active. Then, the dynamical conductivity takes 
a constant (energy independent) value of Gq = e^/(4^), which results in frequency- 
independent absorption of na ~ 2.3 %, where a = /Tic (in cgs units) is the fine 
structure constant. Hence, in this model, the absorption of light in graphene is so- 
called universal, given just by a which in general describes the strength of interaction 
between light and matter. Let us note that the graphene dynamical conductivity 
related to the vertical across the Dirac point transitions is frequency independent, 
a{u)) oc {l/uj)D{uj) = const, simply because the joint density of states D{u) is linear 
with energy (frequency). When considering the linearly polarized light, the contribution 
of individual fc-conserving transitions to the total optical conductivity depends on the 
angle between the given momentum k and the vector of light polarization [51], but the 
universal conductivity is restored after integration over all angles. 

Transforming Gq into measurable quantities, such as the optical transmission T and 
reflection R, one finds that i? <C T (i? = n^a^T/A |16]), and therefore, the absorption 
of light in graphene, and also the fine structure constant a itself, can be estimated in 
a straightforward transmission experiment, T = 1 — na. This fact was confirmed in 
recent experiments in the visible range on the self-standing graphene membrane |49j . 
see Fig. HI and further explored by other groups [60l [61] . The universahty of the optical 
conductivity in graphene has been also discussed in the context of epitaxial graphene 
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Figure 5. (a) Real part of the optical conductivity (per one valley) as a function of 
the Fermi level (in units of an arbitrary energy Eq) calculated using the self-consistent 
Born approximation for short-range scattering centers [50J, reprinted with permission 
of the Journal of the Physical Society of Japan, (b) Schematic diagram of intra- and 
inter-band transitions in doped graphene. The non-vanishing Fermi level gives rise to 
the characteristic onset of optical conductivity at 2Ef- Intra-band transitions results 
in the appearance of the Drude-peak centered at zero frequency. 

on SiC substrate [62] in a wide range of photon energies. 

Nevertheless, we should mention that the declared universality of the graphene 
absorption is only approximative (what unfortunately limits its applications in 
metrology), since Hamiltonian ^ represents just the simplest model of the graphene 
band structure and therefore neglects the non-linearity of bands at high energy distance 
from the Dirac point. In other words, the absorption of light in graphene is as "universal" 
as "ideal" are massless Dirac fermions in this system. A theoretical analysis of these 
corrections was performed by Stauber et al. [16], who discussed not only deviations from 
non-hnearity due to finite width of 7r-bands as given by Eq. ([T]), but also included the 
effects of next-nearest hopping integral 7o, see Fig. [2j Experimentally, the deviations of 
the optical conductivity from the "universal" value Gq? have been probed in the visible 
range by Fei and coworkers IM] . Another departure from "universality" of the graphene 
transmission is provided by saturation effects, due to which graphene can find its use as 
a saturable absorber [63l [Ml ISSl 166] . 

It is worth noticing that a nearly universal conductivity has been also reported by 
Kuzmenko et al. [67] for bulk graphite (conductivity per single atomic sheet). Naively, 
one can seek for the origin of this effect in a relatively weak inter-layer coupling 
in graphite, but a deeper theoretical analysis is also available: the appearance of 
the universal conductivity (again per sheet) in multi-layer graphene stacks has been 
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Figure 6. The experimentally determined (real part of) optical conductivity of 
graphene placed on Si/Si02 substrate for different positions of the Fermi level |69) . 
Reprinted by permission from Macmillan Publishers Ltd: Nature Physics [69] . 
copyright (2008). 

considered by Min and MacDonald [68j in relation to the presence of the chiral symmetry 
in structures with AB and/or ABC stacking sequences. 

A more general theoretical approach to the optical conductivity of graphene, 
which includes the effects of carrier density, temperature, as well as carrier scattering 
(both short- and long-range scatterers considered), has been presented already in 
"pre-graphene" era by Ando et al. |50] in the framework of the self-consistent Born 
approximation, see Fig. [5^. Later on, further calculations were performed by Peres et 
al. IB] as well as Stauber et al. [57]. A simple analytical expression for the optical 
conductivity was found by Gusynin et al. |52] followed by Falkovsky and Varlamov |55] . 
These works account for both inter-band transitions which are blocked at energies 
below twice the Fermi energy {2Ef) as well as intra-band transitions which result in 
the appearance of the Drude peak centered at zero frequency (see Fig. [5]). As can 
be expected, these models uncover the universal conductivity Go = e^/ (4^) in the 
limit of high energies. If the Fermi energy approaches the neutrality point, the optical 
conductivity is universally Go = e^/(4^) at any finite frequency but singular at w = 0. 
The apparent value of a{uj = 0) related to the so-called minimum conductivity is one 
of the lively discussed issue in the graphene physics, without clear consensus between 
different experiments and theoretical models so far [501 El IH 1301 EOl l3T] . 

Experimentally, the optical response of exfoliated graphene flakes placed on Si/Si02 
substrates has been investigated in the mid-infrared range in the reflection configuration 
by Wang et al. [72J who reported the pronounced modification of graphene's optical 
properties as a function of the position of the Fermi level. Further experimental data 
taken in both reflection and transmission mode, accompanied by a detailed analysis 
have been presented by Li et al. [69] and later on also by Mak et al. [60] . The spectrum 
of the optical conductivity extracted from these experiments (see Fig. |6]) is in overall 
agreement with theoretical predictions. However, some of the observed features, such 
as the striking scaling of the onset of inter-band transitions with the gate voltage or the 
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intriguing finite absorption in between the Drude peak and the 2Ep onset remain to be 
firmly clarified. Possibly, these observations might be due to effects of disorder and/or 
of electron-electron interactions, which as yet may not be fully accounted in theoretical 
models [53] . 

3.3. Magneto-spectroscopy 

Cyclotron motion of charge carriers and the related cyclotron resonance phenomenon 
(absorption of light at cyclotron frequency cuc) is primarily a classical effect, probably 
the most representative for the magneto-optical spectroscopy. Such motion is not only 
characteristic for a conventional charged (e) particle with mass m, which precesses with 
the frequency of uc = eB/m. The solution of the classical equation of motion for 
charged particle with energy e, depending linearly on momentum p [e = vpp), also 
results in the cyclotron motion but with the frequency uc = eB /{\e\/vp), in which one 
easily identifies the energy dependent mass m = \e\/vp. This latter equation, equivalent 
to the Einstein relation between mass and energy, invokes again the relativistic-like 
character of electronic states in graphene. Perhaps surprisingly, the strictly speaking 
classical, i.e., linear with the magnetic field, cyclotron resonance has not been clearly 
observed in graphene so far, thought the existing experimental results are well explained 
by quantum mechanical approach. 

In a quantum-mechanical picture, the application of the magnetic field B 
perpendicular to the graphene plane transforms the continuous electronic spectrum into 
discrete and highly degenerate Landau levels (LLs) [73] : 

En = sign{n)vF\/'2\e\fiB\n\ = sign{n)Ei'\J\n\, = 0, ±1, ±2 . . . (3) 

whose positions are defined by a single material parameter, the Fermi velocity vp 
{El = vrV^heB). The de generacy of each Landau level is C{B) = gygs\eB\/h, where 
we take into account both spin and valley degeneracies. This LL spectrum consists of 
electron levels {n > 0), hole levels (n < 0) and the zero LL {n = 0) which is shared 
by both types of carriers and which is responsible for the unusual sequence of the 
quantum Hall effect in graphene [21 H]. We also immediately see that LLs in graphene 
are non-equidistant and they evolve as \/B, see Fig. [7^, what both can be understood 
as a consequence of the extreme non-parabolicity (in fact linearity) of the bands. The 

Table 1. Table of various graphene- related quantities together with their numerical 
values. The energy e and the Fermi level Ep in the last column are expressed in meV 
and the Fermi velocity vp in units of 10^ m.s"""^. niQ denotes the bare electron mass. 



Quantity 


Relations 


Numerical values 




Effective mass 


m = lel/vj, 


1.76 X lQ-^{\e\/vl) 


mo 


Density of states 


gsgv\e\/i2Trvjph'^) 


1.47 X 108(|£|/4) 


cm^^.meV""'^ 


Carrier density 


gsgvEj./iiTTvltt^) 


7.35 X l{f{Ep/vpf 


cm^^ 


Energy of 1^' LL 


El ^ VFy/2heB 


m.ivF^jB[T] 


meV 
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Figure 7. (a) A characteristic \/S-dependence of LLs in graphene shown for a few 
low-index levels (b) Dipole-allowed inter-LL transitions for two different filling factors 
V — 1 and 6. Real parts of the (low-temperature) optical conductivity (proportional 
to the optical absorption) determined using Eq. ^ are plotted in panels (c) and (d) 
for filling factors v = 1 and 6, respectively. The grey curve is calculated for the same 
position of the Fermi level as in ;y = 6 case, but at temperature of 150 K. 

unusual -y/B-dependence of LLs in graphene is responsible for its surprising sensitivity 
to the magnetic field. Experimentally, the well-defined LLs have been spotted in this 
system down to 1 mT and almost up to temperature of liquid nitrogen [71] . It might be 
realistic that Landau level quantization in pure graphene could also be observable in the 
magnetic field of the Earth (-Bsarth ~ 10~^ T), which is unique for a condensed-matter 
system. 

Interaction of light with graphene in a quantising magnetic field has been explored 
several times both theoretically and experimentally over the few past years [HI [751 
[76l|73|78l|73|80l[8llll8lll71|82l|72IH3]. To describe the magneto-optical response 
of graphene, we start with a single-particle model and use the Kubo-Greenwood 
formula [HU |85l [86]. Then, the longitudinal conductivity takes the following form: 

Re(a±(u;,5)) = ^^^El("^l^±Hr7^^ (4) 

where < /„ < 1 stands for the occupation of the n-th LL and 7 is a phenomenological 
broadening parameter. The matrix elements of the velocity operators -0+ = (vx—ivy) / \/2 
and {)_ = (t)a; + iVy)/V2 are \{m\v+\n)\'^ = a6\m\,\n\+i and \{m\v^\n)\'^ = a6\m\,\n\-i, 
respectively, with a = f|,/2 if n or m equals zero, and otherwise a = Vp/A [76ll8T| [79 | l87] . 

We see that graphene exhibits a relatively rich (multi-mode) magneto-optical 
response, where energies of individual resonances which correspond to individual inter- 
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LL transitions scale as \/B, preserving this unique property of the Landau level fan 
chart in graphene. Let us now discuss a few basic facts which imply from Eq. @: All 
dipole- allowed inter-LL transitions in graphene follow the selection rules |?t,| — )■ |r;,| + 1 
and \n\ \n\ — 1, which are active in the and a~ polarization of the incoming light 
[HUlHIj, respectively. The possible transitions can be divided into three groups (j > 1): 
Inter-band resonances L_j — > Lj+i and L^j^i — )■ Lj at energy Ei (^^/J^fT + ^/]^ , intra- 

band resonances Lj — )■ L^+i and L_j_i — )■ L_j with energy Ei (^^/J^fT — y/j^ and the 
mixed L_i(o) ^o{i) resonance, involving the n = LL, having energy of Ei. From 
a purely terminological point of view, only intra-band resonances can be referred to 
as cyclotron resonance (CR), nevertheless, the term Cyclotron resonance experiments 
is used frequently to describe any of inter- or intra-band inter-LL transitions in the 
current literature, and becomes thus equivalent to a more general term of Landau level 
spectroscopy. 

A scheme of optically active transitions for a hypothetical system of 2D Dirac 
fermions at two, arbitrary chosen, filling factors u = 2 and 6 is illustrated in Fig. [7)d. 
The corresponding optical-conductivity spectra are presented in Figs. [7b and [Hi. The 
intra-band transitions appear at low energies and are followed by inter-band resonances 
at higher energies. There is however no distinct separation between these two types of 
transitions, what is in contrast to the case of conventional 2D systems based on gapped 
semiconductors, but somehow similar to the case of narrow-gap II/VI compounds 
structures [HHl [H9]. Nevertheless, in graphene, we deal with only one type of atomic 
orbital, and therefore both intra- and inter-band transitions follow similar selection 
rules: modulus of the LL index is changed by 1. This is again in contrast to the case of 
convectional 2D system, made for instance of GaAs, for which the inter-band transitions 
conserve the LL index and their dipole moment is due to different s- and p-orbitals of 
the conduction and valence band, respectively. Owing to the electron-hole symmetry of 
the graphene band structure, two different inter-band resonances, such as, for example 
L_2 — ^ L3 and L_3 — )■ L2 in Fig. [7b, niay appear at the same energy. Such degenerated 
in energy transitions are however active in opposite circular polarization of light. At low 
temperatures one may expect at most two different intra-band transitions but a series 
of inter-band transitions. The situation is even more complex at higher temperatures, 
when the thermal spreading of the Fermi distribution exceeds the separation between 
Landau levels. The intra-band absorption (CR) may then also reveal a multi-mode 
character, as seen, for example, in Fig. [TJi, where the grey curve depicts the thermal 
activation of Lq — J-Li and L2 — J-La transitions. Such multi-mode intra-band absorption 
spectrum, which envelop corresponds, nota bene, to the classical cyclotron resonance 
discussed at the beginning of this section, was recently observed in graphene (on the 
surface of bulk graphite) by Neugebauer et al. [71]. 

Another interesting point is that the velocity operators (and also their matrix 
elements) are magnetic-field-independent in graphene, contrary to the case of standard 
systems with a parabolic band [90]. In consequence, if the occupation difference {fn — fm) 
between the initial and final state Landau levels is not changed with the magnetic field 




Figure 8. (a) Typical multi-mode magneto-transmission spectrum of graphene, taken 
on a few-layer epitaxial specimen (3-5 layers) prepared on C-terminated surface of 4H- 
SiC [76] . The inset shows a schematic of the evolution of Landau levels with applied 
magnetic field, and possible optical transitions, (b) Fan chart of observed inter-LL 
transitions, showing clear \fB evolution of all absorption lines. Reprinted from [76], 
copyright (2006) by The American Physical Society. 



(which is often the case of inter-band transitions), the oscillator strength of such a 
transition rises as a/B, see Eq. (jl]) with uj ~ a/B. The transitions involving the n = 
Landau level are also interesting in this respect. In the quantum limit, once the Fermi 
energy is pinned to the n = level, the oscillator strength of the degenerate (non- 
polarized) L_i(o) — 7- Lo(i) transition does not depend on the particular occupation /q 
of this level. This is due to the fact that the change in one of the occupation factors 
(/_i — /o) and (/o — /i) relevant for the L_i — )■ Lq and Lq — )■ Li transitions, respectively, 
is compensated by the same gain for the other one. For the same reason, the oscillator 
strength of the unpolarized transition involving the n = LL increases as a ^J~B at high 
magnetic fields. This is very much in contrast to the case of conventional electron gas 
with parabolic dispersion relations, for which the oscillator strength of the cyclotron 
resonance is just a measure of the electron concentration, including the high-field limit. 
However, the doping of graphene in this limit can still be deduced from the comparison 
of the intensities of the (polarized) Lq — )■ Li and L_i — )■ Lq transitions. 

Experiments: Magneto-optical measurements, mainly magneto-transmission, have 
been up to now performed on three different kinds of graphene samples. First 
experiments were carried out on epitaxial graphene prepared by the thermal 
decomposition of the surface of silicon carbide [THl IHD] , see Fig. [HI To be more specific, 
these CR data have been obtained on a multilayer epitaxial graphene (MEG) specimen 
prepared on C-terminated surface of 4H-SiC. The Dirac-like spectrum, genuine of the 
graphene monolayer, found to dominate in MEG structures, has been initially an 
intriguing observation but this is today a well-established experimental fact [291 l9T] - 
The electronic bands in MEG grown on the carbon face of SiC are alike those of a 
single layer because of preferentially rotational and not Bernal-type layer stacking in 
this material [92] l9T] . To our knowledge, no reliable magneto-optical measurements 
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Table 2. Fermi velocity in various graphene specimens deduced by different 
experimental techniques: magneto-optics (MO), ARPES, and STS. Values for nearest- 
neighbour coupling parameter 70 were obtained via relation 70 — 2hvp/ {\/3ao), 
7o[eV] = 3.088wF[10^m.s-i]. 



Graphene specimen 


VF (10^ m.s-i) 


70 (eV) 


Remark 


Exfoliated, Si02 substrate 


1.12 


3.46 


MO 


m 


Exfoliated, Si02 substrate 


1.09 


3.37 


MO 


m 


Epitaxial (C-terminated surface) 


1.02-1.03 


3.15-3.18 


MO 


IZilHl] 


Epitaxial (C-terminated surface) 


1.13 


3.49 


STS 




Epitaxial (C-terminated surface) 


1.0 


3.1 


ARPES [9l] 


Decoupled on surface of bulk graphite 


1.00 


3.09 


MO 


m 


Decoupled on surface of bulk graphite 


0.79 


2.4 


STS 


[23 



have been reported on samples grown on Si-face of 4H-SiC (or 6H-SiC), where mostly 
the standard (well-known from bulk graphite) Bernal stacking is set. On the other 
hand, Si-face samples have been used in a number of ARPES experiments, see 
e.g. |23l EU |25]. The magneto-optical experiments on epitaxial samples were soon 
followed by measurements on a single flake of exfoliated graphene [THl HE] , see Fig. [9l and 
recently, also by CR in decoupled graphene sheets on the surface of bulk graphite [71] . 
To compare, whereas the large-size samples (up to a few mm^) covered with MEG 
flakes (up to a few mm^) allow for the relatively easy and precise CR measurements 
on well-defined, practically undoped graphene sheets (no ~ 5 x 10^ cm~^ [22]), the 
demanding experiments on single flakes of exfoliated graphene are more subtle. They 
have been successfully realized with a differential (modulation) technique when using 
the gate voltage to tune the Fermi energy and therefore to block the apparent transitions 
in the reference signal. The magneto-transmission experiments offer several pieces of 
important information about electronic properties of graphene, we will summarize them 
in the following. 

Landau level fan charts and Fermi velocity: A clear illustration of the characteristic 
^J\Bn\ scaling of Landau levels, in fact equivalent to the observation of linear dispersion 
relations of carriers, has been the first important feat of the Landau level spectroscopy 
of graphene systems [76]. Sadowski et al. [76] observed a practically perfectly '^\Bn\- 
scaled Landau level fan chart in MEG structures and extracted vp = 1.03 x 10^ m.s~^ 
for the Fermi velocity in this system. Jiang et al. [78J and subsequently Deacon et 
al. [H] found somehow higher values vp ~ 1.1 x 10^ m.s^^ for gated (exfoliated) 
graphene flakes on Si/Si02. We note that Landau levels in MEG structures as well 
as in graphene flakes floating on the graphite surface have been recently also visualized 
using tunnelling spectroscopy in magnetic fields [281 129]. The Fermi velocity found 
in these later STS experiments agrees well with magneto-transmission data in case of 
MEG structures [761 US] but it is surprisingly low 0.79 x 10^ m.s~^ [28j in case of 
graphene flakes on graphite, for which the magneto-spectroscopy gives much higher 
value of 1.00 x 10^ m.s~^ [74j. Individual values of the Fermi velocity determined in 




Energy (meV) „(io«cm') 

Figure 9. (a) Magneto-transmission spectra of exfoliated graphene on Si/Si02 
substrate obtained by Jiang et al. [78]. Individual intcr-LL transitions are identified 
in the inset. (b) Density dependence of the two-contact resistive voltage and 
photoconductive response of an exfoliated graphene specimen at the fixed magnetic 
field and fixed energy of CO2 laser, for details see Deacon et al. Reprinted 
from [7H1I1H]! copyright (2007) by The American Physical Society. 



various samples by different experimental techniques have been summarized in Tab. [2j 
Beyond simple band models: Relatively small but noticeably deviations of bands 
from their ideal linearity, of the order of a few percent at large ±0.5 eV distance from 
the Dirac point, have been found by a combination of far and near- infrared magneto- 
optical experiments performed on multi-layer epitaxial graphene by Plochocka et al. ^7] . 
These deviations were revealed by a departure of the observed transitions from a simple 
A/B-dependence, which rises with the photon energy of the probing light. No signs of 
the electron-hole asymmetry has been found in these experiments. On the other hand, 
the traces of the electron-hole asymmetry have been reported by Deacon et al. [IE] 
in exfoliated graphene placed on Si/Si02 substrate, who estimated the difference in 
the electron and hole Fermi velocities to be of the order of few percents. Magneto- 
transmission experiments if carried out on neutral graphene specimens may also bring a 
relevant information on a conceivable appearance of a gap at the Dirac point. Working 
in the limit of low magnetic fields, Orlita et al. have estimated a gap to be smaller 
than 1 meV in quasi-neutral MEG structures [82] and its maximum possible value of a 
fraction of 1 meV in graphene flakes on graphite substrates [74] . 

Scattering/disorder : Cyclotron resonance measurements on graphene, in particular 
in the limit of low magnetic fields (low frequencies) can be effectively used to estimate 
the scattering time and/or mobility of carriers. For instance, Orlita et al. [HS] (working 
in fields down to 10 mT range) have shown (see Fig. [TO]) the possibility to achieve the 
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room-temperature carrier mobility exceeding 250 000 cm^/(V.s), which is the record 
value among all other known materials. Recent predictions based on the analysis of 
transport data taken on exfoliated graphene [93] have been thus confirmed. Moreover, 
the closer look at the CR lineshape can offer valuable information about the scattering 
mechanisms. An increase of the (LL-independent) linewidth nearly following the 
A/S-dependence indicates that the short-range scattering dominates in this particular 
case [91]. Using also the magneto-optical methods, the mobility of charged carriers 
exceeding 10^ cm^/(V.s) up to the temperature of liquid nitrogen has been determined 
for high-quality graphene flakes on the surface of bulk graphite [H]. CR measurements 
thus became a straightforward method to characterize the quality of graphene specimens, 
see, e.g., Ref. [95] . 

Electron- electron interaction: Since the discovery of graphene, the effects of 
electron-electron interaction were a subject of particular interest in this material. 
Nevertheless, a great majority of experimental results obtained on various graphene 
systems are fairly well understood within single particle models. This also concerns a 
number of magneto-transmission studies [761 El ESI [71] . Characteristically, they display 
a regular, defined by a single parameter vp series of transitions which are thus very 
tempting to be assigned as those between single particle Landau levels. However, the 
excitations between highly degenerate Landau levels are known as nontrivial processes 
which involve the effect of electron-electron interaction. Such an electron-hole excitation 
is characterized by its wavevector (nota bene proportional to electron-hole separation). 
The specific shapes of dispersion relations of inter and intra Landau level excitations 
are central for the many-body physics of the integer [961 El II] and fractional [971 IHHl [99] 
quantum Hall effects, respectively. We know from this physics that, when considering 
a single parabolic band of a conventional 2DEG (with equidistant Landau levels), the 
energies of optically active k = inter Landau level excitations correspond to those of 
single particle excitations. This can be viewed as a consequence of Kohn (or Larmor) 
theorem or seen as a result of perfect cancellation of Coulomb binding and a positive 
exchange term for the k = electron-hole excitation. This reasoning does not hold 
for a 2D gas of Dirac electrons, for which the exchange term may even largely exceed 
the Coulomb binding, both in addition expected to be different for different pairs of 
Landau levels. The apparent approximate validity of the Kohn theorem in graphene is 
a surprising effect, and in our opinion calls for further clarifications on the theoretical 
background. The first theoretical works dealing with this problem have been already 
published [100[ \101\ 1102] . We note however here that although small but noticeable 
deviations from a perfect single particle scaling of inter Landau level excitations have 
been already reported in experiments on exfoliated graphene structures [78]. Very 
recently, Henriksen et al. [83] have reported changes in the energy of the L_i(o) — )■ Lo(i) 
transition, especially pronounced at high magnetic fields, when tuning the Fermi energy 
in between n = —1 and n = 1 Landau levels. Both these observations [78] 183] are 
discussed in terms of electron-electron interactions but perhaps also include some effects 
of disorder. Notably, magneto-optics allows for probing the nature of quasi neutral 
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Figure 10. Part (a): Far infrared transmission T plotted as — In T as a funetion of 
the magnetic field at T = 2.0 K. The dashed lines denote the expected transitions for 
vp = 1.02 X 10^ m.s~^. The inset shows the transmission spectrum at i? = 0.3 T. 
Part (b): FIR transmission measured at T = 2 K and low magnetic fields. Successive 
spectra are shifted vertically by 0.05. The part (c) shows the peak position, width and 
area for the L-i(o) ^Lq(i) transition. The dashed line in part (c) is a least squares fit 
to the peak positions. Reprinted from ^82^, copyright (2008) by The American Physical 
Society. 



graphene in high magnetic fields (a possible appearance of a gap in the zero LL at filling 
factor zero), which recently became a subject of many theoretical considerations and 
experimental works [35[ 1103] . 

4. Bilayer graphene 
4.I. Band structure 

The interest in the physics of a graphene bilayer started after magneto-transport 
experiments of Novoselov et al. [20] who showed a characteristic quantum Hall effect and 
a Berry's phase of 27r with magneto-resistance measurements of this material. Carriers 
in a graphene bilayer are characterized by a finite mass but are clearly distinct from 
those in conventional systems with parabolic bands and referred to as massive Dirac 
fermions. Similarly to the case of graphene, the basics of the electronic band structure of 
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Figure 11. Crystal structure of the Bernal-stacked graphene bilayer with the 
corresponding S WM hopping parameters 70 , 7i , 73 and 74 . 



the graphene bilayer have been estabhshed in 1947 by Wallace [5] and further developed 
in the band structure model of bulk graphite introduced by Slonczewski, Weiss and 
McClure HOSl [106] (SWM model). As a matter of fact, it is not graphene, but 
the graphene bilayer, which is the basic unit in the construction of Bernal-stacked bulk 
graphite. 

The presence of four atoms in the unit cell of the graphene bilayer implies the 
appearance of four bands in the vicinity of the Fermi level, instead of two in the 
graphene monolayer. The SWM model, reduced to a true bilayer, represents a solid 
basis for their description. This model implies one intra-layer 70 and three inter-layer 
7i;73;74 hopplug integrals, see Fig. [TTl In addition, the difference in on-site energies of 
stacked and unstacked triangular sublattices A', as well as, the bias voltage U applied 
across the layers are usually taken into account when modelling the band structure of 
the graphene bilayer. 

As only 70, 71 and U parameters predominantly affect the shape of bands, we neglect 
the others for a while and express the Hamiltonian of the graphene bilayer as [2T] : 

/ U/2 vpn^ \ 

fr_ -U/2 VF-K , 

VF-K^ -U/2 71 ' 

V vpi: 71 U/2 I 

which corresponds to the case of atoms in layers 1 and 2 ordered as (Ai, i?2, -Bi). 
The electronic bands obtained by the diagonalization of this Hamiltonian have simple 
form of: 

-E'1,2 = "-£'4,3 = 

= - (71V2 + f/V4 + vIp' ± ^|l^/^ + vlp\-^l + u^)) . (6) 

These bands are symmetric with respect to the zero energy level, defined by the touching 
point of the bands E2 and £"3 (so-called charged neutrality point) for U = 0, see Fig. [T2k. 
A nearly parabolic shape of bands E2 and E^ close to the charge neutrality point 
allows us to introduce an effective mass m = '-fi/{2vp) with a relatively low value of 
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Figure 12. A schematic view of the band structure in graphene bilayer, showing its 
modification via bias voltage U apphed across the two layers, (a) J7 = and (b) U ^ 0. 
The arrows show dipole-allowed transitions in this system which determine its optical 
response in the infrared spectral range. 

m ^ O.OSmo [107] comparable to a mass found in narrow gap semiconductors. Assuming 
no doping or external gating, the Fermi level is located just at the charge neutrality 
point and similarly to graphene, the graphene bilayer can be viewed as a zero- gap 
semiconductor (semimetal). As sketched in Fig[T2b. the external bias U applied across 
the layers, transforms this gap-less band structure into a system characterised by a finite 
energy gap E^^p ^ U. 

4-2. Optical conductivity at B = 

The relevant optical response of the graphene bilayer is expected in the infrared 
range, see the illustrative Fig. [121 ^is the characteristic band energy scale is given 
by the parameter 71 ^ 0.4 eV [13]. Experiments have followed the well established 
theoretical background [871 I109[ I110[ llll[ 1112] . Two groups almost simultaneously 
reported on experiments performed on a single gated flakes of the exfoliated graphene 
bilayer |113[ I108[ I114[ 1115] . see Fig. [131 overall confirming the expected electronic band 
structure given by Eq. ([6]). A more closer data analysis revealed a weak electron- 
hole asymmetry, going beyond the Hamiltonian ([5]) and quantified by parameters 
74 ^ 150 meV and A ^ 18 meV |113] . see also Tab. [3l No pronounced influence 
of the trigonal warping, expressed predominantly by the 73 parameter, was found in 
these experiments. Very recently, analogous experiments have been performed on more 
sophisticated devices with two gates incorporated, allowing independent tuning of the 
carrier density and voltage U [116[ I117j . A possibility to smoothly tune the energy 
gap from zero up to E'gap ~ 200 meV has thus been demonstrated. The possibility 
of relatively easily tuning of the gap in the graphene bilayer depicts the interesting 
potential of this system in designing the electronic devices. 
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Figure 13. Real part of the optical conductivity of the bilayer graphene as a function 
of the gate voltages, i.e. at various hole (a) and electron (b) densities measured by 
Li et al. |108) . The inset in part (a) schematically shows the band structure of the 
graphene bilayer without and with parameters 74 and A' taken into account. Position 
of individual band-edge inter-band transitions in the optical conductivity is shown in 
the inset of the part (b). Reprinted from [108] . copyright (2009) by The American 
Physical Society. 



4.3. Magneto-spectroscopy 

In most of the experimentally studied cases, the graphene bilayer evokes the quantum 
mechanical character of its band structure under the magnetic field applied across the 
2D plane. Following the Hamiltonian ([5]), the energy ladder of Landau levels with indices 
n > 0, is given by: 

rpn rpn 

-^1,2 ~ -^4,3 ~ 

= - (71V2 + in + l/2)El ± ^^t/A + (n + 1/2)^^? + Ef/A) (7) 

The eightfold {8eB/h) degeneracy of the LL with n = is twice the degeneracy of 
all other levels and this results in a characteristic quantum Hall effect with the Berry 
phase of 27i [SOj [21]. At energies around the neutrality point, the Landau levels are 
formed of E2 and E^ bands, and their energies are practically linear with B. As shown 
in Fig. [T^ . a departure from (saturation of) this linear dependence appears however 
at higher energy distances from the neutrality point (from the zero LL at energy zero) 
and/or for Landau levels with higher indices. 

Similarly to the case of the graphene monolayer, the magneto-optical response of 
graphene bilayer can be described within the Kubo-Greenwood formalism, employed 
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Figure 14. a) LLs in graphene bilayer, evolving nearly linearly with B at lower 
fields, but bending to sub-linear dependence at higher B and/or higher energies, (b) 
Dipole-allowed inter-LL transitions for the case of two filling factors v = A and 8. 



Table 3. Band structure parameters of the graphene bilayer 71, 74 and A' derived from 
infrared optical experiments performed on the exfoliated graphene bilayer specimen. 



Band structure parameter 


Value 


References 


71 


400 meV 


|113lll()8j 




378 meV 






« 350 meV 




74 


150 meV 


|ll3lll08j 




140 meV 


[niins] 


A' 


18 meV 


|113lll08j 




22 meV 


|114llll5j 



Table 4. Various quantities related to the graphene bilayer in the parabolic 
approximation of E2 and bands. The numerical values are calculated for 
71=375 meV and vp = 10^ m.s~^. toq denotes the bare electron mass. 



Quantity Relations Numerical values 

Effective mass TO ^i/{2v\) 0.033 toq 

Density of states 5sfft,7i/(4w|.7r?i^) 2.76 x 10i° cm'^.meV"^ 

Carrier density 5^5^71 |£'F|/(4w|,7r?i2) 2.1ii ^ IQ^^\Ef\ cut'^ 

Cyclotron energy TiWc 2heBvp/ji 3.51i?[T] meV 
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in the previous Chapter, but the individual matrix elements are not easy to evaluate. 
Up to now, this problem was approached approximatively by Mucha-Kruczyhski et 
al. |118] and also some numerical results are available |119[ 1120] . In general, we expect 
that all transitions between LLs arising from all four bands are dipole-allowed, if the 
selection rules n — )■ n ± 1 are obeyed, and they are active in the cr^ polarization of 
the incoming light, respectively. A few transitions in the vicinity of the Fermi level, 
including their polarization, are schematically shown in Fig. [Hb, for two arbitrarily 
chosen filling factors f = 4 and 8. Certain additional complications arise from the 
appearance of the trigonal warping, expressed mainly by the parameter 73, which is 
not included in the Hamiltonian ([5]) and which induces the second set of dipole-allowed 
transitions between LLs. Such transitions appearing between levels whose indices differ 
by 3iV ± 1, where = 1, 2, 3 . . . [i87j, may have their oscillator strengths comparable 
with n — )■ n ± 1 transitions but only at low magnetic fields ^87j . 

The scheme of inter Landau level transitions in a graphene bilayer is largely 
simplified when a parabolic approximation of E2 and bands is considered [21], i.e. 
when the Hamiltonian ([5]) is reduced to 2x2 matrix (the Ei and E^ split-of bands are 
neglected). This approach seems to be fairly sufficient when considering the low-energy 
excitations in an undoped or weakly doped graphene bilayer. The LL spectrum then 
takes the form of E'^ = —E2 = hcOc\Jn{n + 1) for ri > 0, or equivalently: 

En = sign(n)/iu;c\/|^|(|^| + 1), n = 0, ±1, ±2 . . . , (8) 

where stands for the cyclotron frequency Uc = eB/m with m = 7i/(2t>|,). In this 
approximation, the Landau levels are perfectly linear with the applied field B. Let us 
note that in the limit of high \n\ (practically even for \n\ > 1), the LL spectrum Eq. ([8]) 
has the form En ~ sign(n)/ia;c(|^| + 1/2), typical of conventional massive particles. 

Within the parabolic approximation, the matrix elements of the velocity operators, 
|(m|-0±|n)p, to be set in the Kubo-Greenwood formula, have been evaluated by Abergel 
and Fal'ko [87]. They obtained: Kml-O^ln)!^ = Ib'-^c + l) <^|m|,|n|±i (Jot m 7^ 

and n 7^ 0), where Ib = \J'h/ (eB) denotes the magnetic length. Transitions involving 
doubly degenerated n = LL level must be taken into separately. The spectrum is 
obviously expected to depend on the actual position of the Fermi energy and special 
care has to be taken to account for possible splitting of the n = Landau level (including 
the elucidation of the character of this splitting). To prevent confusion, we note that 
our convention of LL indexing differs from the one used by Abergel and Fal'ko, whose 
Landau level spectrum reads: E^ = —E2 = TujJc\Jn{n — 1) for n > 1, as compared to 
our Eq. (E]). 

The available magneto-optical data on a graphene bilayer are not as rich as on 
graphene. Probably the only report is by Henriksen et al. |107j . who examined a 
single gated flake of an exfoliated graphene bilayer in (differential) far infrared magneto- 
transmission experiments. In this work [107] the authors investigated the transitions 
between the adjacent LLs (intra-band transitions) as a function of the magnetic fleld 
and the carrier density (fllling factor). A relatively good agreement with the simplifled 




Figure 15. Magneto-transmission spectrum of an exfoliated graphene bilayer 
measured by Henriksen et al. |107j . Upper left inset: the quantum Hall effect of the 
bilayer graphene sample measured in situ at = 18 T. Upper right inset: schematic of 
the zero- field dispersion relation of bilayer graphene. Reprinted from [107] , copyright 
(2008) by The American Physical Society. 



LL energy spectrum ([71) has been obtained for low carrier densities. At higher filhng 
factors, deviations appear due to the finite energy gap induced by the gate voltage, as 
recently explained by Mucha-Kruczyhski et al. |121] . 

5. Graphite 

5.1. Band structure 

The renewed interest in the properties of bulk graphite is a direct consequence of 
the outbreak of the current graphene physics. As a 3D crystal, graphite is a system 
characterised by a higher degree of complexity compared to graphene, nevertheless, 
both material share many common properties. 

The physics of bulk graphite has been reviewed several times, see e.g. Refs. |l2l US] . 
Here we mostly focus on optical properties which distinctly uncover the massless and 
massive Dirac-like electronic states also in this material. The appealing possibility to 
trace the "relativistic" carriers not only in graphene monolayer and bilayer, but also 
in bulk graphite (which is definitely easier to handle!) resulted in a number of works 
which offer new pieces of information, new interpretations of old data but often also 
rediscoveries of well-established knowledge. 

Starting with the pioneering work of Wallace [5], the fundamentals of the 
graphite band structure have been formulated by Slonczewski, Weiss and McClure in 
fifties |105[ I104[ 1106] , as already mentioned in Section HI This model describes the 
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Figure 16. Brillouin zone of graphite and its electronic band structure along the 
H-K-H line. Whereas electrons in the vicinity of the K point behave as massive Dirac 
fermions with mass twice enhanced in comparison to the true graphene bilayer, the 
holes around the H point have nearly linear dispersion and behave as (nearly) massless 
Dirac fermions in graphene, but with an additional double degeneracy. 




Figure 17. Segment of the crystal structure of A-B stacked graphite with coupling 
constants used in the standard SWM model. 



band structure along the H-K-H line of the 3D Brillouin zone, see Fig. [161 which 
is responsible for most of the electrical and optical properties of bulk graphite. The 
SWM model mostly implies six hopping integrals 70, ... , 75, as sketched in Fig. [TTl and 
additional parameter A, usually referred to as pseudogap and related to the difference 
of the crystal field on atom-sites A and B. By definition, this pseudogap is related to a 
similar parameter in an isolated graphene bilayer as A' = A — 72 + 75. 
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To describe optical response of graphite, we leave the full complexity of the 
SWM model and follow a simplified approach taking account only two most relevant 
hopping integrals 70 and 71, which describe intra- and interlayer tunnelling, respectively 
[122[ 1119] . In this way, the band structure of bulk graphite along the H-K-H is obtained 
by the diagonalization of the Hamiltonian: 
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(9) 



which is formally equivalent to that of an unbiased graphene bilayer ([5]), with the 
same ordering of the atomic wavefunctions, but with an effective coupling tuned by the 
momentum in the direction perpendicular to layers, A = 2cos(7rA;2). The appearance 
of the cosine band in kz, having the amplitude of 271, thus reflects the periodic 
crystal ordering along the c-axis, similarly to formation of minibands in semiconductor 
superlattices |123j . 

The straightforward diagonalization of the Hamiltonian ([9]) gives the band structure 
in the form of Eq. ([6]), i.e. electronic bands equivalent to the graphene bilayer, but with 
an effective coupling A71: 

Ei,2 = -E,,s = - ((A71) V2 + vIp' ± v^(A7i)V4 + t^V(A7i)') • (10) 

Hence, depending on their momentum along the c-axis, electrons in bulk graphite behave 
as massive Dirac fermions in graphene bilayer, but with an effective coupling A71, which 
directly gives their effective mass m = A71/ {2vp). For example, electrons at the K point 
[kz = 0) behave as massive Dirac fermions with an effective mass enhanced twice (A = 2) 
in comparison to true graphene bilayer. Carriers (in real graphite holes) at the H point 
{kz = 0.5) have a character of massless Dirac fermions due to the effectively vanishing 
inter- layer coupling (A = 0). Bulk graphite thus shares some common properties with 
the graphene monolayer as well as its bilayer. 

5.2. Optical conductivity at B = 

In reflection and/or transmission experiments at zero magnetic field the whole optical 
response corresponds to an average over all momenta kz along the H-K-H line. With 
this averaging, the concluded response (per one sheet) is surprisingly close to the optical 
conductivity of an isolated graphene monolayer, as pointed out by Kuzmenko et al. [67] , 
see also Refs. |125] 1126] . Hence, a nearly universal optical conductivity can also be 
observed in bulk graphite. 

5.3. Magneto-spectroscopy 

Historically, the verification of the band structure expected in the framework of the 
SWM model and estimates of individual tight-binding parameters have been often done 
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Figure 18. Positions of inter-LL transitions related to the H point of bulk graphite 
observed by Toy et al. [124] in magneto-reflection experiment, see Fig. [TBI The observed 
well-defined v^-dependence is characteristic of massless Dirac fermions. The splitting 
of the line with the lowest energy is induced by existence of the pseudogap A. Reprinted 
from [124) . copyright (1977) by The American Physical Society. 



using optical experiments in magnetic fields. In this case, the main contribution to the 
optical response is provided just by the K and H points, where Landau bands become 
flat, what leads to singularities in the joint density of states |119l 1120] . Cyclotron 
resonance of massive electrons around the K point was reported in several works, see 
older papers Refs. |127[ I128[ I129[ 1130] and also recent measurements [125] , showing well- 
defined response which scales nearly linearly with B. Much less information has been 
collected about the massless Dirac fermions around the H point but nevertheless the 
features in magneto-reflection spectra which follow the -y/S-dependence, typical of 2D 
massless particles, have been already reported by Toy et al. |124j . This work |124] may be 
considered as probably a very first observation of massless Dirac fermions in solids. Later 
on, magneto-transmission experiments on a very thin graphite specimen (fslOO sheets) 
allowed for a deeper analysis of the H point optical response |13Hll32j . Magneto-optical 
response of the H point of graphite includes a set of transitions equivalent to those of 
graphene [76l[78| H8]. but it is nevertheless significantly richer |133] . Since there are four 
carbon atoms in the unit cell of graphite (and only two in the unit cell of graphene) the 
degeneracy of Landau levels in graphite (at the H point only) is doubled as compared 
to the case of graphene. This results in a set of additional dipole allowed inter Landau 
transitions at the H point of graphite which are indeed observed in experiments, see 
transitions marked with small Greek letters in Fig. [191 

The basic understanding of the complex magneto-optical response of bulk 
graphite |134j . see Fig. [191 is possible just by viewing this material as an effective 
graphene monolayer and an effective bilayer with a coupling strength enhanced twice 
in comparison to a true graphene bilayer. Interestingly, this model including only 
two material parameters, 70 and 71, is capable to catch the basic (magneto-)optical 
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(a) Energy/B"=(meV/T"=) (b) B"^(T"^) (c) B(T) 



Figure 19. (a) Transmission spectra of a thin graphite layer at selected magnetic 
fields. The plotted energy is scaled as E/\/B to emphasize the massless Dirac 
fermion-like features originating at the H point (indicated by dashed vertical lines). 
Arrows denote transitions arising at the K point which evolve (nearly) linearly with 
B. (b) Positions of the absorption lines related to the H point as a function of 
\fB. The solid and dashed lines represent expected positions of absorption lines for 
vp — 1-02 X 10^ m/s (c) Positions of the absorption lines related to the K point 
as a function of B. The solid lines show the expected dipole- allowed transitions 
in a graphene bilayer with an effective coupling 271 calculated using Eq. (1101) for 
VF = 1-02 X 10^ m/s and 71 — 375 meV. The inset schematically shows the observed 
inter-band transitions in the effective bilayer. Reprinted from jl34) . copyright (2009) 
by The American Physical Society. 

properties of bulk graphite, a complex 3D material. This creates an interesting link 
between physics of bulk graphite and graphene monolayer and bilayer. The values of 
band structure parameters determined by magneto-transmission experiments |134] are 
7o = (3.20 ± 0.06) eV and 71 = (375 ± 10) meV. The presented model is however 
evidently simplified, and therefore, let us now consider the role of other tight-binding 
parameters 72 ■ ■ ■ 75 and A. 

The integral 72, due to the hopping between next-nearest graphene planes (see 
Fig. [T7|) . is dominantly responsible for the presence of electron and hole pockets around 
K and H points and therefore, i.e., for the semimetallic character of graphite. This is 
clearly evidenced for instance by magneto-transport experiments |135l 11361 11371 1138] . 
Using the language of the effective graphene bilayer and monolayer, we obtain an 
effective n-doped bilayer at the K point with the Fermi level Ep ~ 20 meV and similarly, 
p-doped monolayer around the H point with Ep ^ —20 meV. It is important to note 
that magneto-optics and magneto-transport probes are sensitive to very different states 
selected from those along the K-H dispersion line. Magneto-optics is selective with 
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respect to the states at K and H points, because it is sensitive to the maxima in joint 
density of states (integrated along the direction) of Landau bands which becomes 
flat at the extreme of the Brillouin zone along the z direction. In turns, magneto- 
oscillations in transport or de Haas-van Alphen experiments result from the subsequent 
coincidences between the Fermi energy position and the maximum in the density of 
states. Those latter maxima appear not only at K and H point but importantly also at 
the intermediate point somewhere in the middle of the dispersive line. It appears that 
K- and intermediate-point maxima, both reflecting rather massive (bilayer) character, 
are responsible for the two characteristic sequences of magneto-oscillations observed in 
magneto-transport experiments on graphite |139l 11381 1140] . 

Trigonal warping of the band structure of graphite is described in the SWM model 
by the 73 parameter (73 ^ 300 meV |l3]), and its relative influence on the band 
structure is scaled with the parameter A |120] . Therefore, the band structure is mainly 
modified in the vicinity of the K point and its influence at the H point is completely 
cancelled. Probably, the most pronounced effect related to trigonal warping is the 
observation of the harmonics of the basic CR mode due to the iT-point electrons. 
Interestingly, 73 parameter negligibly affects the LL energies at the K point, but it 
induces new selection rules allowing transitions between LLs differing by 3A^ ± 1 in their 
indices |127[ I14H I129[ 1130] . whose relative strength is increasing with the decreasing 
magnetic field [87]. Let us note, however, that determination of LLs within the full 
SWM model is only possible via direct numerical diagonalization of a truncated infinite 
matrix |142] . 

The effects of 74 and A parameters are analogous to the the case of a true bilayer, 
i.e., they in principle induce the electron- hole asymmetry of the band structure as 
discussed in the previous Section. Again, the effective strength of 74 scales with A, 
giving the most pronounced effects at the K point, where a relatively weak asymmetry 
has been observed |143l 11441 [GTj 1134] 1145] . At the H point, the parameter A corresponds 
to a (pseudo)gap opened at the Dirac point in the graphene-like band structure and is 
estimated in the range of several meV |124l 11311 H5l 1132] . 

6. Few-layer Bernal-stacked segments of graphite 

A fairly established knowledge about optical properties of the 2D graphene monolayer 
and bilayer, as well as 3D bulk graphite logically leads to a question about the properties 
of graphene stacks composed of three and more layers with Bernal ordering of sheets. 
As pointed out by Partoens and Peeters )146] the 10-layer graphene stack resembles 
already graphite very much, but the investigation of the graphene-to-graphite evolution 
is interesting and instructive. To get insight into the electronic states of a few-layer 
graphite stack, we use a relatively simple and intuitive model. We start with the 
electronic band structure of bulk graphite and select discrete values of momentum kz 
along the c-axis, = (7r/2))[l — m/{N + 1)], with m = 0,2,...(A^ — 1) for odd 
or with m = 1,3,... (A^ — 1) for even A^, for which a standing wave appears in a 
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Figure 20. The expected dynamic conductivity of Bernal-stacked A^-graphene layers, 
as calculated within a simple band structure model of these systems discussed in the 
text. All graphene stacks are assumed to be neutral. The dynamic conductivity of 
the monolayer is constant and equals ire^/h. The bilayer conductivity is calculated 
according to the theory developed in Ref. [87], and assuming 71 = 385 meV. The 
conductivity of the higher order stack is the sum of the conductivities of its effective 
bilayers and if applicable of a monolayer (if N is odd). For example: the trilayer 
conductivity is the sum of the conductivity of the monolayer and the effective bilayer 
with the coupling parameter of I.4I71; the conductivity of the = 6 stack is the sum 
of the conductivities of three effective bilayers with effective coupling parameters of 
0.4571, 1.2571, and I.8O71. 



stack with layers. Using a mathematical language, this approach corresponds to 
a decomposition of the Hamiltonian (assumption that 72 = 75 = is compulsory) 
describing a graphite stack with A^-layers (dimension 2N x 2A^) into a set of sub- 
Hamiltonians (4 x 4). For even N, each segment corresponds to the Hamiltonian of 
the bilayer with the effective coupling constant Am7i, which is directly related to the 
allowed discrete values of momentum as = 2cos/€m,. For odd number of stacked 
layers A^, we get very similar results, (A^ — l)/2 effective bilayers and one monolayer-like 
Hamiltonian (2x2). Hence, in this way, the problem of the band structure of few-layer 
graphite is reduced to the already discussed band structures of graphene monolayer 
and bilayer, by simple introducing the appropriate effective coupling parameters. This 
approach has been discussed by Partoens and Peeters |122j and further developed by 
Koshino and Ando [1471 11191 1120] who also included the influence of the magnetic field. 
Perhaps more elaborated, but less intuitive approach can be found, e.g., in Refs. |148[ H5]. 

Experimentally, only very first measurements which probe the optical properties of 
graphene stacks with N > 2 have just appeared for stacks up to A^ = 8. They display 
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a clear evolution of the absorption spectra as a function of and well resemble the 
features which can be seen in the traces of the calculated optical conductivity, expected 
in these structures (see Fig. 1201) . 

7. Summary and Outlook 

Optical spectroscopy has in recent years been successfully applied to investigate the 
properties of new two dimensional allotropes of carbon and to revise the properties of 
bulk graphite. Majority of these experiments have been devoted to establish or verify 
the band structure character of the investigated materials and they certainly helped to 
establish and/or to better understand the characteristic nature of massless or massive 
Dirac electronic states in graphene [691 EHl SHI SHI [60], in graphene- like MEG layers on 
SiC [761 ig |82] , in bilayer graphene [IIHl HOHl [113 [iH] , or in graphite [ISH |67l [134] 
and its few layer segments [149] . Characteristically for these semimetalic materials, the 
methods of infrared spectroscopy have been frequently used but also (nearly) visible 
optics appeared as a relevant tool to uncover their properties [IHl [HI [EU I150j . The 
authors of this review are particularly attached to magneto-spectroscopy techniques 
and consider them among the most effective methods to studying the band structure 
parameters. Estimations of the Fermi velocity in graphene based systems or of the 
band curvature (effective mass) in bilayer graphene structures are the primary results 
of applications of these methods, see Tabs. [THU In complement to convectional electric 
transport measurements, the magneto-optics is also effective to evaluate the efficiency of 
carrier scattering [821 [95] and becomes a unique method if the electric contacts cannot 
be applied. 

More challenging for experimentalists, such as the authors of this review, is to 
uncover the unexpected properties. Identification, with magneto-optical experiments, 
of graphene-like bands in MEG structures [76] was certainly a surprise, but is now a 
well acknowledged fact [921 [291 [91] . Similarly, the surprisingly extremely weak scattering 
efficiency anticipated for graphene from magneto-optical experiments founds also its 
justification in a very recent theoretical work [151] . The deviations from simple models, 
tentatively attributed to the effects of electron-electron interactions, observed in the 
magneto-optical studies of graphene |78l [83] or in studies at i? = ^69l [60] are equally 
appealing. 

It is obviously risky to foresee the further developments in optics of graphene- 
related structures but probably safe to say that they will be very much dependent on 
the progress in sample preparation: fabrication of easy handled large area samples, 
of systems with higher electronic quality and certainly new structures. Optical and 
magneto-optical methods will continue to be used to extract the relevant band structure 
parameters of these materials. These methods are likely to be of particular importance 
with respect to intense efforts of band gap engineering in different graphene systems [20] 
(e.g., graphene bilayer under the apphed voltage [1161 11171 11211 1115] or chemically 
functionalised graphene [1521 11531 1154] ). The problem of the specific, for linear 



Dirac electronic states in graphene systems: 



Optical spectroscopy studies 



31 



bands, role of electron-electron interactions in inter Landau level transitions is to be 
clarified both on the experimental and theoretical ground. Better quality samples 
could certainly help in experiments, as they recently did in the case of the first 
observations of the fractional quantum Hall effect in graphene [981 EH]- Intriguingly, 
the classical collective phenomena such as low energy plasmons are absent in the 
experiments so far. This is to be verified, for example, in low-field low-frequency 
magneto-optical experiments on samples with well defined shape, as already done for 
standard 2D semiconductor systems |155[ 1156] . Low-energy plasma excitations might 
be very unique in graphene |157l 11581 11591 I160[ 1161] and also interesting from the 
viewpoint of THz applications |162] . The predicted nonlinear effects related to cyclotron 
resonance and/or plasmon excitations are also to be examined |163[ 1164] . A possibility 
to detect the cyclotron resonance emission which may perhaps easily appear in the 
system with non-equidistant Landau levels, due expectable inefficient inter-Landau level 
Auger processes |165[ I166[ 1167] , is worth to verify. This could give rise to efficient {B 
tunable) sources of infrared radiation or even to THz lasers |168] . (Magneto-)optical 
experiments may be also decisive in confirmation of the Dirac subbands [169] when 
fabricating the artificial graphene |170j (by honeycomb lithography of a conventional 
2DEG). Approaching the illusive quantum electrodynamics effect of Zitterbewegung is 
perhaps another challenge for optics of graphene [iTl 11711 11721 1173] . 

It should be stressed that this review covers only the basic linear transmis- 
sion/reflection optical studies of graphene-based structures. For example, it does not 
include the vast domain of Raman scattering experiments |174l 1175] . We just briefly 
note, that Raman scattering measurements on graphene systems are so far limited to in- 
vestigations of the phonon response. However, since we deal with resonant processes, the 
phonon response also includes information on the electronic bands and, for example, its 
analysis has become one of the most popular method of primary characterization of dif- 
ferent graphene structures. The electron-phonon interaction |176[ 1177] is another impor- 
tant issue that has been successfully studied in graphene using Raman scattering tech- 
niques |178l 1179] , including also the application of magnetic fields |180l 11811 1182] . Obser- 
vation of Raman scattering signal due to electronic excitations in graphene systems |183] 
would be highly interesting in view of investigations of single particle versus collective 
plasmon excitations and/or dispersion relations (due to electron-electron interaction) 
of inter Landau level excitations. Another class of optical experiments, which are not 
discussed here, concerns time resolved studies and the related investigations of carrier dy- 
namics. There are already few reports in this area [182 [Ml HHSl IMl [Wl [1881 [Ml IHO] , 
which is likely to expand in the future. Investigations of carrier dynamics with time 
resolved studies in the far infrared range and under application of a magnetic field might 
be interesting with respect to search for far infrared emission from graphene structures. 
One may finally speculate that these are perhaps the unexpected discoveries, such as 
the recent report on broad band emission from graphene in the visible range [150] (in- 
duced by oxygen plasma etching), which will open the entirely new domain in optics of 
two-dimensional allotropes of carbon. 
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